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Abstract: By exploiting the analyticity and boundary value properties of the thermal Green 
functions that result from the KMS condition in both time and energy complex variables, we treat 
the general (non-perturbative) problem of recovering the thermal functions at real times from 
the corresponding functions at imaginary times, introduced as primary objects in the Matsubara 
formalism. The key property on which we rely is the fact that the Fourier transforms of the 
retarded and advanced functions in the energy variable have to be the "unique Carlsonian analytic 
interpolations" of the Fourier coefficients of the imaginary-time correlator, the latter being taken 
at the discrete Matsubara imaginary energies, respectively in the upper and lower half-planes. 
Starting from the Fourier coefficients regarded as "data set" , we then develop a method based on 
the Pollaczck polynomials for constructing explicitly their analytic interpolations. 



1. Introduction 

In the standard imaginary-time formalism of quantum statistical mechanics (tracing back to Mat- 
subara |l5|] ) and, later on, of quantum field theory at finite temperature (see e.g. jl4| and references 
therein) , there arises the a-priori non-trivial problem of recovering the "physical" correlations at 
real times starting from data at imaginary times. More specifically, the correlations at imaginary- 
time of observables (or, more generally, of boson or fermion fields) in a thermal equilibrium state at 
temperature T — /3 _1 are defined as periodic (or antiperiodic) functions of period f3, and therefore 
they are equivalently characterized by their discrete mode expansion 4 Q n exp(— iC, n T) in terms 
of the so-called "Matsubara energies" iC n , where ( n — 2mr/P (or (2n + l)ir/ (3). 

The problem of recovering the correlations at real time, or equivalently the retarded and ad- 
vanced Green's functions at real energies, from the previous sequence of Fourier coefficients {G n } 
admits a unique and well-defined theoretical solution in terms of the notion of "Carlsonian analytic 
interpolation of this sequence" . This can be achieved || , and will be recalled below in Sect. [| if the 
imaginary-time formalism is embedded in the (conceptually more satisfactory) general description 
of quantum thermal states as KMS states jnj. However, as suggested by the lattice approach of 
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the imaginary-time formalism, it may be interesting to have a concrete procedure for constructing 
satisfactory approximate solutions of this problem when one starts from incomplete data sets. 

In this paper we give a precise algorithm for the previous reconstruction problem; this mathe- 
matical method is presented in Sect. Moreover, in the subsequent Sect. |J, the method is applied 
to the case when the data are finite in number and affected by noise. 

Let us consider the algebra A generated by the observables of a quantum system. Denoting 
by A,B,... arbitrary elements of A and by A — ► A(t) (A = A(Q)) the action of the (time- 
evolution) group of automorphisms on this algebra, we now recall the KMS analytic structure of 
two-point correlation functions (A(ti)B(t2)) q , in a thermal equilibrium state Dp of the system 

at temperature T = ft ■ 

By time-translation invariance, these quantities only depend on t — t± — t% , and we shall put 

W AB {t) = {A{t)B) Qp , (1) 
W' AB (t) = (BA(t)) Q)3 . (2) 

In finite volume approximations, the time-evolution is represented by a unitary group e lH so 
that 

A{t)=e iH>t Ae- iH '\ (3) 

where H' = H — fiN, H being the Hamiltonian, fi the chemical potential, and N the particle 
number; under general conditions, the operators e~@ H have finite traces for all ft > (see e.g. JTTf). 
Then the correlation functions are given, correspondingly, by the formulae 

W AB (t) = j-Ti {e-^' A(t) B} , (4) 
W' AB {t) = ^-Tr{e-^'BA(t)}, (5) 

where Z p = Tr e _/3H ' . 

One then introduces the following holomorphic functions of the complex time variable t + ij: 

Gab (t + iy) = ^Tr \ e -^ H ' A(t) e^ H ' b) , (6) 
Zp l J 

analytic in the strip {t + ij; t G K, —ft < 7 < 0}, and 

G' AB (t + iy) = j-Tr {e-^ H 'Be-^'A(t)} , (7) 

analytic in the strip {t + 27; t € R, < 7 < ft}, which are such that: 

lim G AB (t + i 1 ) = W AB (t), (8) 

7— >o 

lhnG' AB (t + i^)=W AB (t). (9) 

7>0 

From (0) , (Q) and the cyclic property of Tr , we then obtain the KMS relation 

W AB (t) = Tre- f3H 'A{t)B = Tr B e' 011 ' A{t) = G' AB {t + ift), (10) 
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which implies the identity of holomorphic functions (in the strip < 7 < (3) 

G AB (t + i( 7 -p)) = G' AB (t + i 7 ). (11) 

According to the analysis of |ll[] in the Quantum Mechanical framework and of || in the Field- 
theoretical framework, this KMS analytic structure is preserved by the thermodynamic limit under 
rather general conditions. 

In the case when the algebra A is generated by smeared-out bosonic or fermionic field operators 
(field theory at finite temperature), the principle of relativistic causality of the theory implies 
additional relations for the corresponding pairs of analytic functions (G, G"). In fact, this principle 
of relativistic causality is expressed by the commutativity (resp. anticommutativity) relations for 
the boson field <£(x) (resp. fermion field S'(x)) at space-like separation: 

[${t,x),${t',x')] =0 (resp. {&(t, x),V(t',x')} = 0) for (t - t'f < (x - x'f . (12) 

In this field-theoretical case, we can choose as suitable operators A the "smeared-out field op- 
erators" of the form A = f$(yo 7 y)f(yo,y)dy dy (resp. / &(y ,y)f{y ,y)dy dy), where / is 
any smooth test-function with (arbitrary small) compact support around the origin in space- 
time variables. For the observable B, we can then choose any operator A x obtained from A by 
the action of the space-translation group (which amounts to replace the test-function f(yo,y) 
by fivoiV) = f(ya,y — x ))- It then follows from ( |l2] ) that the corresponding analytic functions 
GAAvit + ij) and G' AA (t + ry) (satisfying (pi])) have real boundary values Waa^) and (t) 
which satisfy, on some interval \t\ < t(x, /), coincidence relations of the following form: 

WAA„(t) = W AAw (t) in the boson case, (13) 
WAA„(t) = -W'aa ( f ) in the fermion case. (14) 



Then, in view of identity (|l 1\) , the coincidence relations ([L3|) and ([14|) imply the existence of a 
single analytic function Qaa^ (t + 17) which is such that: 



a) in the boson case: 



b) in the fermion case: 



Qaa. = Gaa„ for - 13 < 7 < 0, (15) 
Saa. = G' AA ^ for < 7 < 0; (16) 



Qaa„ = G AA „ for - /3 < 7 < 0, (17) 
Gaa. = -G'aa„ for < 7 < /3. (18) 

Correspondingly, it follows that Qaa^ is either periodic or antiperiodic with period if3 in the full 
complex plane minus periodic cuts along the half-lines {t + 17; t > t(x,f), 7 = k/3, k € Z} and 
{t + ir, t < -t(x, /), 7 = k/3, k € Z}. 

These analytic functions Qaa^ (t + *7) are smeared-out forms (corresponding to various test- 
functions /) of the thermal two-point function of the fields <P (or in the complex time variable. 
In other words, this thermal two-point function can be fully characterized in terms of an ana- 
lytic function Q(t + 17, x) (with regular dependence in the space variables) enjoying the following 
properties: 

a) Q(t + ij,x) — eQ(t + 1(7 — (3), x), where e = + for a boson field, and e = — for a fermion field; 

b) for each x, the domain of Q in the complex variable t is C \ {t + 17; \t\ > |cc|; 7 = k(3, k £ Z}; 
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c) the boundary values of Q at real times are the thermal correlations of the held, namely: 

lira Q(t + ij,x) =W{t,x), (19) 

7— *o 

hm Q(t + ij,x) = W'(t,x), (20) 

7>0 

where in hnite volume regions, W and W" can be formally expressed as follows (a rigorous 
justihcation of the trace-operator formalism in the appropriate Hilbert space being given in ^|): 

W{t, x) = \- Tr e-? H <P(t, x) <P(0, 0), (21) 

Z P 

W'(t,x) = — Tie' (3H <P(0,O)<P(t,x), (22) 
Zp 

for the boson case, and similarly in terms of ^(t, x) for the fermion case. 
In this analytic structure, we shall distinguish two quantities that play an important role: 

i) the restriction G{i~i, x) of the function Q to the imaginary axis is a /3-periodic (or antiperiodic) 
function of 7 which must be identified with the "time- ordered product at imaginary time^\ 
considered in the Matsubara approach of imaginary-time formalism. In the latter, this quantity 
or its set of Fourier coefficients plays the role of initial data. 

ii) The "retarded" and "advanced" two-point functions 

R(t,x) = i6(t)[W(t,x) -eW'(t,x)l (23) 
A(t,x) = -i6(-t)[W(t,x) - eW'{t,x)}, (24) 

which are respectively the "jumps" of the function Q across the real cuts {t; t > \x\} and 
{t; t < — \x\}. These kernels have an important causal interpretation; in particular, R describes 
the "response of the system" to small perturbations of the equilibrium state. The knowledge of 
R and A and, consequently, of W — W = — i (R — A) allows one to reconstruct W and W by the 
application of the Bose-Einstein factor 1/(1 — e T,3w ) to their Fourier transforms W{uj) 
(this procedure being an implementation of the KMS property in the energy variable lu). 

The rest of the paper is devoted to the problem of recovering the "real-time quantities " R and A, 
starting from the "time-ordered product at imaginary times" as initial data. This will require the 
conjoint use of the analytic structure of Q in complex time and of its Fourier-Laplace transform 
in the complex energy variable. In fact, the key property on which our reconstruction of real-time 
quantities relies is the following one: the Fourier-Laplace transforms R and A of the functions 
R and A, which are defined and analytic respectively in the upper and lower half-planes of the 
energy variable u>, are analytic interpolations of the set of Fourier coefficients {Gn} of the function 
Q at imaginary times, the latter being taken at the Matsubara energies u> — i( n . Moreover, the 
uniqueness of this interpolation is ensured by global bounds on R and A, according to a standard 
theorem by Carlson ||. The basic equalities that relate R(iCn) an d A(i£„) to the corresponding 
coefficients Gn will be called "Froissart-Gribov-type equalities" for the following historical reason. 
A general n-dimensional mathematical study of the type of double-analytic structure encountered 
here has been performed in [[| in connection with the theory of complex angular momentum, 
where the original Froissart-Gribov equalities had been first discovered (in the old framework of 
S-matrix theory). The fact that this structure is relevant (in its simplest one-dimensional form) 
in the analysis of thermal quantum states has been already presented in |^| in the framework of 
Quantum Field Theory at finite temperature. 
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2. Double Analytic Structure of the Thermal Green Function and 
Froissart Gribov type Equalities 

In the following mathematical study we replace the complex time variable t + ij of the introduction 
by t = i(t + 17) in such a way that, in our "reconstruction problem" treated in Sects. |^ and ||, the 
initial data of the function G(t, •) considered below correspond to real values of t. Up to this change 
of notation, this general analytic function Q(t, •) can play the role of the previously described two- 
point function of a boson or fermion field at fixed x. However, since the only variables involved in 
the forthcoming study are r and its Fourier-conjugate variable C, the extra "spectator variables" , 
denoted by the point (•), may as well represent a fixed momentum (after Fourier transformation 
with respect to the space variables) or the action on a test-function / (as for the correlations of 
field observables A — A(f) described in the introduction). 

Let us summarize the analytic structure that we want to study. 

Hypotheses. The function G(t, ■), (r = U + iv, u,V G K), satisfies the following properties: 

a) it is analytic in the open strips k[3 < u < (k + l)/3 (v € K, k € Z, (3 = 1/T) and continuous 
at the boundaries; 

b) it is periodic (antiperiodic) for bosons (fermions) with period [3, i.e. 

Q( T+ R .) = |^( r ' ') for bosons, (tGC), 
^ ' 1 ' \ — Q(t, ■) for fermions, (t€C); 

c) sup \G(u + iv,-)\ < C\v\ a , (v el; C, a constants). (25) 

-k/3<u<(k+l)(3 

We shall treat both the boson and fermion field cases at the same time by exploiting the 2/3- 
periodicity of the function Q{r, •). To this purpose, we take the Fourier series (in the sense of 
L 2 [-/3, [3]) of G(t, ■), which we write 

-. +00 

G ^^ = W E Qn{-)e-<"\ C« = ~n, (26) 

n— — 00 

and whose Fourier coefficients are given by 

Qn(-)= / 9{T,-)e^ T dT. (27) 
J-p 

It is convenient to split expansion (^6|) into two terms as follows: 

-. +00 

S {+) (r,-)^-J25i +) (-)e- 1 ^, g(+)(-) = a„(-), (n = 0,l,2,...), (28) 

ZP n=0 

^(^•) = ^ £ ^-K-)e-^, OtK)=enU, (n=-l,-2,...), (29) 



then. 



Ql +) {-)= f Q {+ \r,-)e^dT, (n = 0,1,2,...), (30) 
J-p 

Q { n\-)=l G ( - ] {r,-)e^dT, (n = -l,-2,...). (31) 
J-P 
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We now introduce in the complex plane of the variable r = u + it) G R) the following domains: 
the half-planes X± = {r £ C |Imr ^ 0}; the "cut-domain" I + \ £+, where the cuts S+ are 
given by S+ = {r £ C | r = k(3 + iv, v > 0, k £ Z}, and 21 \ SI, where S- = {r £ C | r = 

o 

+ ro, i? < 0, k £ Z}. Moreover, we denote by ^4 any subset A of C which is invariant under the 
translation by k(3, k £ Z (e.g. S± s 2±\ S±, etc.) (see Ref. Accordingly, the periodic cut- 

r-plane C \ (.5*+ U S-) will be denoted by 7J T . We now introduce the jump functions Jy^pAv, •) 
and J(kp\(v, •) that represent the discontinuities of Q^ + \t, •) and G^~\t, •) across the cuts located 
respectively at Rer = u = k(3, v > 0, and at Rer = u = k/3, v < 0, (k £ Z): 

J[k0)( v r) = +ilKn{g(+\ki3 + e + iv,-) - g(+\kj3 - e + iv,-)\ , (v > 0, fc G Z), (32) 
J^ } (v,-) = -i\^\G^ ) (kf3 + e + iv,-)-G { ~ ) (k[3-e + iv 1 -)\, (» < 0, fc G Z). (33) 

6>0 

Let us note that these definitions are well-posed and appropriate because, as we shall see in the 
following theorem, C?W(r, •) and Q^ir, ■) are holomorphic in the cut-domains Z_ U £"+] and 
X + U [Z_\ S-]j respectively. Moreover, we suppose hereafter that the slow-growth condition fl25| ) 
extends to the discontinuities J^\(v, ■), that turn out to be "tempered functions" Q. Finally, in 
view of the periodicity properties of G(t, •), it is sufficient to consider only the strip, in the r-plane, 
defined by -a < u < 2(3 - a (0 < a < (3), v £ R (see Fig. |). 

We then introduce the Laplace transforms of the jump functions across the cuts located at 
Rer = 0, and at Rer = (3; i.e. 

r+oo 

j[o) (C, ') = J o J$(v, e~ (v dv, (C = £ + iv, ReC > 0), (34) 



J ( ( o ) ) (C,-)= J { (0) \v,-)e~^dv, (ReC<0), (35) 

J — oo 

/■+0O 

J$(C, •) = / •) e-t° dv, (ReC > 0), (36) 



^(C ') = / 4? K ') e ~ Cw ( Re C < 0). (37) 



J —CO 

We can state the following theorem. 

Theorem 1. If the functions Q (r ; ■) andJ^^Av, ■) satisfy the slow-growth condition uniformly 

o o o 

iti 77 T = C \ [S+ U ) Up i/ie closure, the following properties hold true: 

i) The function Q^ + \t, •) (respectively Q^-~\t, •)) is holomorphic in the cut-domain X_ U £*+] 

o 

(respectively X+ U £"-])• 

ii-a) T/ie Laplace transforms jft (Cj*) an< ^ (d') are holomorphic in the half-plane ReC, > 0. 
TTie Laplace transforms J^g?((, ■) and J^m(Ci ■) o^e holomorphic in the half-plane ReC, < 0. 
ii-6) J^(£, •) and J(j3){Ci ') belong to the Hardy space H 2 MC|^\ where C = G C | i?e£ > 

<5, <5 > e > 0}. ■7(0)'' (Cj") an< ^ J(p)(G>') belong to the Hardy space H 2 ^C^j, where = {C G 

C |iteC<M>e>0}. 
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Fig. 1. Integration paths used in the proof of Theorem 



iii-a) In the case of the boson statistics the symmetric combinations Cj' + ' b )(£, ■) — J(n)(Ci') + 



(Ci ') and '^(C)') — J(c>) + J(8'\{Gi') interpolate uniquely the Fourier coefficients 



.= ?(")/ 



?(-) 



(0) 



'03) 



'(0) 



'08) 



§2„ (•) and Q^m (') respectively (hereafter the superscript (b) stands for the boson statistics). Let 
Cm — 2mir/(3, then the following Froissart-Gribov-type equalities hold: 



(Cm,") 



()) (Cm, ') + J(p) (Cm, ') 



(Cm? ') + J(p){Cm-> ') 

(-), 



-r.(-) 



(■), (m= 1,2,3,...), 
S&>(-), (m = -1,-2,-3,. 



(38) 
(39) 



iii-b) In the case of the fermion statistics the antisymmetric combinations G^ + ^\Q, ■) := J(n\(Ci ') 



J(ff\(C)') an d '^(Cj ') •= ^fo^(C'0 — J(8)iC>') "interpolate uniquely the Fourier coefficients 



.- ?(-)/ 



r(-) 



(0) 



'(« 



'(0) 



'08) 



^2m+i(') an d ^2m+i(') respectively (hereafter the superscript (/) stands for the fermion statistics) 



Let £ m = (2m + 1)~k / (3, then the following Froissart-Gribov-type equalities hold: 

g^f\Cm, ■) = J^(Cm, ■) - 4 + ?(Cm, ■) = 0&L(O, (m - 0, 1,2,3, . . .), 



(0) V^m; 
- f(-) 



^(Cm, ') — <^(q) (Cm; ') J(g) (Cm, ') 



'(/?) 
(-)/ 



03) 



^2m+l('); 



(m= -1,-2,-3,...). 



(40) 
(41) 
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Proof, (i) In view of the Riemann-Lebesgue theorem, and since £ (+) (t, •) G Li[-0,0\, the Fourier 

coefficients Q^\-) tend to zero as n — > oo. From expansion ( pgj ) we have for all r = u + iv, with 
v < 0: 



, +00 



< 



#Xy»», (42) 



n>0 



where if = |£(r, -)| dr. The series X^n>o e converges uniformly in any domain compactly 
contained in the half-plane Imr < 0. In view of the Weierstrass theorem on the uniformly conver- 
gent series of analytic functions, we can conclude that Q^ + \t, •) is holomorphic in the half-plane 
Imr < 0. By using analogous arguments we can prove that Q(~^(t, •) is holomorphic in the half- 
plane Imr > 0. Furthermore, we know from Hypothesis a) that Q(t, ■) = Q^(r, •) + Q^~\t, ■) is 
holomorphic in the strips k/3 < u < (k + l)/3 (fc G Z, v G R), and continuous at the boundaries of 

the strips. We can conclude that G^ + '(t, •) is holomorphic in the cut-domain X_ U [2+\ S+], and 

Q^~'[t, •) is holomorphic in the cut-domain X+ U [I-\ £-]. 

(ii) Property (ii, a) follows easily from the assumption of "temperateness" of the jump functions |Q . 
For what concerns property (ii, b) we limit ourselves to prove that •) belongs to the Hardy 

space H 2 (C^^, since the remaining part of the statement can be proved analogously. To this 
purpose, we rewrite the Laplace transform (|34|) in the following form: 



+00 



^\vr)e- Sv )e-^dv:=J^ s) (C',-), (ReC>0), (43) 



where Re £' = Re £ — S (5 > e > 0). In view of the slow-growth property of jffl (v, •), we can then 

say that the function jtis (v, •) exp(— Sv) belongs to the intersection L 1 [0, +00) n L 2 [0, +00). Then, 

thanks to the Paley- Wiener theorem, we can conclude (returning to the variable C) that J^m(Ci ') 

belongs to the Hardy space H 2 f(C^^ (see Ref. fl2]|). Accordingly, J^(£, •) tends uniformly to 

zero as £ tends to infinity inside any fixed half-plane Re£ > 5' > 6. In particular, J^(Cn, ■), with 
Cn = nit / (3 (n = 1, 2, . . .), tends to zero as n — > 00. 

(iii) We introduce the integral 1^ defined as follows (this method has been introduced by Bros 
and Buchholz g , and will be developed in a more detailed form in [Q within the general framework 
ofQFT): 

4 +) (C,')= f G {+) (T,-)e^dr, (44) 



where the path 7 encloses both the cuts located at u = 0, v > and at u — (3, v > (see Fig. 
In view of the slow-growth condition ( ^5| ) , this integral is well-defined. By choosing as integration 
path a pair of contours (70, 7/3) enclosing respectively the cuts at u = 0, v > and at u = fl, v > 0, 
and then flattening them (in a folded way) onto the cuts (see Fig. Q), we obtain: 

r+00 r + co 

7 (7oU V )(C.-) = / J$(v,-)e-< v dv + e*<> / J$(vr)e-^dv = J$((,.) + eWj$((m 
' Jo 
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Next, we choose the path r y ac , whose support is: ] — a + zoo, —a] U [—a, — e] U [7^] U [e, /3 — e] U 

[7e ] U [/3 + e, 2/3 — a] U [2/3 — a, 2/3 — a + zoo[, where 7e°' ) and 7^ are half-circles turning around 
the points r = and t — (3, respectively (see Fig. [I]). By taking into account the 2/3-periodicity of 
0(+>(t, ■), we get, for ( = ( n = rnr/P, (n = 1,2, . ..): 

/2/3-a 
g (+) (v)e<^T^(+)(.). (46) 
-a 

Then, from the Cauchy distortion argument, we have Z^Ju-yn (Cn> •) = lim e ^o (Cm ')> that is 

^(Cn, + e*"" •) = Gi +) (-)- (47) 

We now distinguish two cases: 
1) n even: n — 2m, ( m = 2m7r//3 (m = 1,2, . . .); then from (|47|) we obtain equalities ( j38|) 



2) n odd: n = 2m + 1, Cm = (2m + l)7r//3 (m = 0, 1, 2, . . . ); then from (|47|) we obtain equalities 

10I). 



We have thus obtained two combinations (symmetric and antisymmetric, respectively) that inter- 
polate the Fourier coefficients Gn~\-)- The uniqueness of the interpolation is guaranteed by the 
Carlson theorem [|| that can be applied since jffliG, ■) and (£, •) belong to the Hardy space 

H 2 ^c|^'^. Proceeding with analogous arguments applied to Q^(t,-) equalities j39| ) and (^|) 
are obtained. □ 

In conclusion, we can say that the thermal Green functions present a double analytic structure 
involving the analyticity properties in the r = u+iv and £ = £+irj planes. The 2/3-periodic function 

Q^ + \t, •) (resp. Q^~\t, •)) is analytic in the cut-domain X_ U [2+\ S+] (resp. T + U S-]) m , its 
Fourier coefficients can be uniquely interpolated (in the sense of the Carlson theorem), and are the 
restriction to the appropriate Matsubara energies of a function <?( +,b ~^(C, ■) (resp. G^' b ^^(Ci ■))> 
analytic in the half-plane ReC > (resp. Re£ < 0). It is straightforward to verify that the 

jump function ji^(v, •) coincides with the retarded Green function, and J^\v, •) coincides with 

the advanced one; analogously, putting z£ = u>, we can identify G {+ ' b ' f) (C, •) and <? (-l6-/) (C, ■) 
respectively with the retarded and advanced Green functions in the energy variable u conjugate 
to the real time t. 

3. Representation of the Jump Function in Terms of an Infinite Set of Fourier 
Coefficients 

First let us consider a system of bosons; since n is even, i.e. n = 2m, £ m = (2rmr)//3, (m = 
0, 1, 2, . . . ), we have: 

G {+ - b) fe,-) =2^W(v)e i? f ^dr. (48) 
Next, recalling that G^(t, •) is /3-periodic, we can write also the following Fourier expansion: 

^\r,-) = ^±G^(^)e-^, (49) 

^ m=0 * ^ ' 

(^)=jy^(r,.)e^ d r^ l -G^ •) . (50) 
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Finally, putting (3 = 2ir, formulae (|49[) , ( |5C| ) can be rewritten in the more convenient form: 

S (+) (r,-) = ^-f:Ql^(m,-)e- im ^ (51) 



m=0 
2n 



0(&)W) = / a (+) (T,-)e imT dr = |a (+ ' 6) (m,0 



(52) 



Recalling once again the /3-periodicity of the function Q^ + \t, •), we write now the Froissart-Gribov 
equalities ( |38| ) as 

C?(+' 6 )(m, •) = J ( ( +< V, •) + Jj+f (m, ■) = 2J ( ( +< b) (m, •) = 2C7 ( +J( TO , •), (m = 1, 2, 3, . . . ).(53) 



It is now convenient to introduce an auxiliary function ji b \v, •), defined as follows 

4 b \v,-) = e- v jj+' b \v,-), («GR+), 



'(0) 



and the corresponding Laplace transform: 

r+oo 



J* ( ) (C,-)= / Ji ) (u ! -)e^dw, (C = ^ + i»?,ReC>-H-*,*>e>0). 



(54) 



(55) 



It is straightforward to prove, via the Paley- Wiener theorem, that J* 6 ^(C, ■) belongs to the Hardy 

! [tl+s) = {C G C | ReC > -1 + 5, <5 > e > 0}. Next, the Froissart- 
■n as 

(56) 



space 



2 (^Ul 



where 



Gribov equalities (53) can be rewritten as 

ji h \m,-) = G^{m + l,-) , (m = 0,1,2,...)- 
Then we can prove the following lemma. 

Lemma 1. The function J* '(— l/2 + i?7, •), (77 G K) can 6e represented by the following series, that 
converges in the sense of the L 1 -norm: 



i=0 



ipe(r]) denoting the Pollaczek functions defined by 



(57) 



(58) 



r being the Euler gamma function, and Pi the Pollaczek polynomials ^,\Id^J. The coefficients a 
are given by: 



C( 



m=0 



—i m + 



(59) 
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Proof. The Pollaczek polynomials P^ a \r]), (77 £ R), are orthogonal in L 2 {— oo,+oo) with weight 
function (see Refs. UM): 



(60) 



For a = 1/2, the orthogonality property reads: 



w(ti)P} 1/a ' J (Tt)P?'''(ri)dri = h,v , I w{ v ) = - 



,(1/2) 



(61) 



(in the following, when a — 1/2, we omit the index a in the notation). Next, we introduce the 
following functions, that will be called Pollaczek functions (of index a = 1/2): 



r 



if) p i(v), 



which form a complete basis in L' 2 (— oo, +oo) [13|]. Since J* (Cr) belongs to the Hardy space 
H 2 ^i +l5 ^ , then Ji^ (— 1/2 + ■) (r/ S K) belongs to L 2 (— oo, +oo). Therefore, we may expand 
Ji b) (-l/2 + irj, ') m terms of Pollaczek functions as follows: 



(62) 



rib) 



(63) 



i=0 



where the series at the r.h.s. of ( |63[ ) converges to J* (— 1/2 + irj, •) in the sense of the L 2 -norm. 
From d63|) we get 



1 



+ DC 



(b) 



1 



(64) 



The integral at the r.h.s. of (|64|) can be evaluated by the contour integration method along the 
path shown in Fig. |2[ and taking into account the asymptotic behaviour of the gamma function 
given by the Stirling formula. We obtain: 



m=0 



— i m + 



(65) 



Finally, from (j5q), (|63| ) and (65) the proof of the lemma follows. □ 
From (|55|), when f = — 1/2 + *r/ (77 € M), we have: 



jib) 



■5 +*»?.■) = 



(66) 



The r.h.s. of ( |6"6| ) is the Fourier transform of J* (u, -)e v l 2 . Noting that J* (— 1/2 + 177, •) belongs 
to L 2 {— 00, +00), but not necessarily to i 1 (~ 00, +00), the inversion of the Fourier transform (pq) 
holds only as a limit in the mean order two, and can be written as follows: 



Ji b \v,-)e v ' 2 = l.i.m. , 

?7o— »+oo V 2tt J _ 



jib) 



IT], 



dri) , (we K + ). 



(67) 



Then, we can prove the following lemma. 
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Lemma 2. The function J*(v, )e v l 2 can be represented by the following expansion that converges 
in the sense of the L 2 -norm: 



where the coefficients aj> are given by: 



_ 00 (_-t \m 

m=0 



-i I m + - 



Pi being the Pollaczek polynomials, and the functions <p£(v) are given by 

^i e V2L e (2e- v )e- e ^e- v/2 , 

Li being the Laguerre polynomials. 
Proof. Let us observe that 
1 



r 



irj 



-t t (i n -i/2) dt 



(68) 



(69) 



(70) 



"e-"/ 2 e- irru dv=^U- e V e~ v/2 Y (71) 



where T denotes the Fourier integral operator. Let us note that the function exp(— e~ v )e~ v ' 2 
belongs to S°°(M.), i.e. the Schwartz space of the C°°(R) functions that, together with all their 
derivatives, tend to zero, f or k >| tending to +00, faster than any negative power of \v\. Therefore, 
we can write (see formula (p3)): 



■■Tip, 



dv 



-v/2 



(72) 
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Substituting in expansion (p7j) to the Pollaczek functions their representation w4), we obtain: 



jf) ( -I +i ,,. ) = g c ,{J=^[p«(- i; i)( e -.-- e -./ 2 )]} 



(73) 



Let us now apply the operator T to the r.h.s. of (73). If we exchange the integral operator T 
with the sum, and this is legitimate within the L 2 -norm convergence, we obtain: 



^p>{j=A p 'H)^" /2 )}}=P>{7- 



Pi -i 



dv 



-v/2 



74) 



Finally, recalling formula (|67|), we obtain the following expansion for the function ji b \v, -)e v l 2 : 



e=o v 



i -i 



dv 



(75) 



whose convergence is in the sense of the L 2 -norm. It can be easily verified that Q 

V2P e (e~ e ^'e-^ 2 ) = i e V2L e (2e~ v ) e^V"/ 2 , (76) 

where Lg denotes the Laguerre polynomials. 

It can be checked that the polynomials Ciiv) = i e \[2 Lg (2e~ v ) are a set of polynomials or- 
thonormal on the real line with weight function w(v) = exp(— v) exp(— 2e~ v ), and, consequently, 
the set of functions <^(v), defined by formula (|70|), forms an orthonormal basis in L 2 {— oo, +oo). 

Finally, from @ we obtain: 



oo 

2 u/2 ji fc) (iv) = , £ai{i e y/2L t (2e- v ), 



-v/2 



oo 

} = ^M>^(v), (we 



1=0 



(77) 



where at = ct/\/2~K, and the functions <&g,(v) are given by formula (|70|). □ 

We now introduce the weighted L 2 -space Lf w) [0, +oo), whose norm is defined by: 



l/li 



w(v) \f(v)\ dv 



o 



1/2 



(78) 



w(v) being a weight function which will be specified in the following. Then we can prove the 
following result. 

Theorem 2. The jump function jfct' (v, •) can be represented by the following expansion: 



jl+- b \v,-) = e v / 2 J2^Mv), (»ei + ), 



(79) 



which converges in the sense of the L^^O, +oo) -norm, with weight function w(v) = e (d€ R + ). 
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Proof. We can write: 

L 

jff l, (v)-e' /s £a<*/(») 



£=0 



12 ,[0,+oo) 



2 \ 1/2 

dv 



1=0 



1/2 



rfc 



In view of Lemma |^ we can thus state that: 



lim 

L — >oo 



0, 



(80) 



(81) 



L= ra) [0,+oo) 



that proves the statement. □ 



Consider now a system of fermions. In this case the function Q^ + \t, ■) is antiperiodic with 
period /3. Then, if we put £ m = (2m + l)ir//3 (m = 0, 1,2, . . .) and /3 = 27r, we have the following 
expansion: 



oo 



(+,/) 
) 



?(+./) 

'(2 



m=0 
/.2tt 



m+~, 



a -i(m+l/2)T 



;f (m + •) = jT 0< + > (r, •) e*"+Va>T dr = \ &+>f) ( m + \ ■) 



(82) 



(83) 



Recalling once again the antiperiodicity of §W(t, •), we write the Froissart-Gribov equalities ( jltf ) 
in the following form: 



;(+./) 



m - 



(+,/) 
(o) 



2' 



J, 



(+,/) 

(27T) 



2" 



2 J 



( + :/) 

(0) 



25 



(+,/) 

(2tt) 



(m = 0,1,2, 



(84) 



We can now proceed in a way strictly analogous to that followed in the case of bosons. We put: 
J { J\v, •) = e- v j\+ J) {v, •) and, accordingly, J { J\C, ■) = / + °° J [ J\v, ^e^dv (f = £ + i?7,ReC = 
£ > —1 + 5, <5 > e > 0). Then, the Froissart-Gribov equalities (§3) now read: 



JP [m+-,- \ = g 



?(+,/) 

'(27T) 



(m = 0,l,2,...). 



(85) 



We can now state the following theorem. 



Theorem 3. i) The function J* ' (if},-), (r] € R) can be represented by the following series, that 
converges in the sense of the L 2 -norm: 



(86) 
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where ipi{ij) are the Pollaczek functions defined by formula (pq), and the coefficients dt arc given 
by: 



(87) 



Pi denoting the Pollaczek polynomials. 

ii) The function J* (v, •) can be represented by the following expansion that converges in the 
sense of L? -norm: 

oo 

ji f) (v,-)=J2 b tMv), (veR+), (88) 

1=0 

where the coefficients bj> are given by bi = d^j \p2/K , and the functions $t(v) are defined by formula 
©■ 

Hi) The function {v, •) can be represented by the following expansion: 

oo 

J Jo) fl ("--)=e"E^W' («€M+), (89) 

1=0 

that converges in the sense of the L? w J0, +oo)-norm with weight function w(v) = e~ 2v , (v € R + ). 

Proof. The proof runs exactly as in the case of the boson statistics, with the only remarkable 
difference that we use the Froissart-Gribov equalities ( ]85| ) instead of (|5^) . □ 

X^? (if), •), which is much more interesting from the physical viewpoint. In order to recover 



We can reconstruct, by the use of this method, the function J* (if), •) but not the function 

(0) 

the function J^'^'(ir), •) we must introduce a more restrictive assumption, requiring the function 

^(o) (£' ") = Ja^ ^(o) ( v >') e ~^ V ^ v to k° holomorphic in the half-plane ReC > —7 (7 > 0). 
Accordingly, in place of the temperateness condition ( p5| ) we assume that J^,'^\v, •) belongs 

to I/ 1 [0, +00) D i 2 [0, +00). Here, for the sake of simplicity, we treat only the case of fermions; 
analogous considerations hold true also in the case of the boson statistics. We can thus suppose 
that the singularities of Xi' (C,-), corresponding to the excited states, all lie in the half-plane 
Re£ < —7, 7 being the smallest damping factor of the spectrum (see Refs. jj],[l(]]). If this is the 
case, jtis'^(irj, •) is analytic, and, moreover, belongs also to L 2 (—oo, +00). We can thus state the 
following result. 

Theorem 4. Let us assume that.j^ f \(,-) is a function holomorphic in the half-plane ReC > —7 
(7 > 0); then Xi'^(«7, •) can be represented by the following expansion that converges in the sense 
of the I? -norm: 

00 

^o)' /) (^-)=E d ^W' ( 9 °) 
1=0 
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where ipe(r]) are the Pollaczek functions defined by formula (pq), and the coefficients d' e are given 
by: 



m— 

denoting the Pollaczek polynomials 



(91) 



Proof. The proof is strictly analogous to the one followed for proving equality ( pT\ ) , and successively 
adapted to the fermion statistics in order to obtain expansion (|86[ ) . The only remarkable difference 
is that now in the expression of the coefficients d' e we have the terms Gr^P ( m + \i ') instead of 
G(2„P { m + fi therefore all the coefficients corresponding to m = 0, 1, 2, . . . , are involved in the 
determination of the function j^'^'fir), •). □ 

~/_ j\ 

Analogous methods and results can be worked out for the function JLy (if), •), assuming that 
J} Q , (C, •) is holomorphic in the half-plane Re£ < 7 (7 > 0). We are then able to reconstruct the 
difference J^.'^\ir), ■) — jL-,'^ (ir), ■) which leads to the determination of the "spectral density" Jl7[ . 



4. Reconstruction of the Jump Function in Terms of a Finite Number of Fourier 
Coefficients 

Up to now we have assumed that all the Fourier coefficients are known, and, in addition, that they 
are noiseless; but this assumption is clearly unrealistic. We now suppose that only a finite number of 
coefficients are known within a certain degree of approximation. We focus our attention on the case 
of the boson statistics, and specifically on the results contained in Lemmas |l| and |^, and Theorem |2| 
The case of the fermion statistics can be treated similarly. We can simplify the notation, without 
ambiguity, by putting: 9$$>(m + l, •) = g m , e v ' 2 J?\v, ■) = F*(v), and J^ } ' b \v, ■) - F(v). Then, 



we denote by gm the Fourier coefficients S}^)'^ + •"■> ') wriCn they are perturbed by noise. We 
now assume that only (N + 1) Fourier coefficients are known within an approximation error of 

order e: i.e. |^ - g m \ < e (m = 0, 1, 2, . . . , N). 
We consider the following finite sums: 



N 



m=0 



i m + 



1 



(92) 



Accordingly, we have a 



(0,00) 



see (|69D). We can then prove the following lemma. 



Lemma 3. The following statements hold true: 

2 



(0,00) 



= F, 



*llL 2 [0,oo) 



ii) 



N— >oo 



-oo. 



— C, (C = constant). 



(£ = 0,1,2,...)- 



(93) 

(94) 
(95) 
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iv) If ko(e, N) is defined as 



k (e,N) = maxik G N : E K | 2 - C 



i.e. it is the largest integer such that X^=o 



2 



< C, then 



lim ko(e, N) = +oo. 

iV— >oc 



v) The sum 



(96) 



(97) 



k 



E| 

£=0 



satisfies the following properties: 

a) it increases for increasing values of k; 

b) the following relationships hold true: 



(k G N), 



fe^oo (AH)2 



(2fc) 2iv , {N fixed). 



(98) 



(99) 



Proof, (i) Equality (|9^) follows from the Parseval theorem applied to expansion (|6^) , and recalling 
that F*(v) belongs to L 2 (— oo, +oo). 

(ii) Let us rewrite the sums a^ e,JV ^ as follows: 



JV 



,(e,JV) 



E & ^ 



z m + 



(100) 



where b$ = y/2{-l) m g$ /ml Now, we can write the following inequality: 

N 

„ £ 6 W P/ 



1 



- / | m + - 



> 



ZZzlbttP e [-i(m + ±)] 



b$P t [-i(N + %] 



(101) 



Let us now recall that in the Appendix of Ref. M the asymptotic behaviour of the Pollaczek 
polynomials P(\—i{m + 1/2)] for large values of I (at fixed m) is proved to be: 



Pt. 



-i [ m + - 



(-l)V 



(2£) m . 



(102) 



Therefore, we have: 





2-<m=0 


b$P e [-i (m + i)] 




b$Pt[-i(N+h)] 




U N 







AT-1 

E 



N 



—M) 

ml 



m-N 



(E03) 
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From (pi]), (|l0|) and dTo|) it follows that for I sufficiently large: 



'N 



(21) 



N 



(104) 



Therefore, lim 



= +oo, and statement (ii) follows. 



(iii) We can write the difference af ,0 °^ — a£ as follows: 



i-iy 



(9m-9$)Pt 



+ E 



rn—0 

(-iy 



i m + 



m=N+l 



i \ m + 



(105) 



In view of the fact that the serie s y^2 X^m=o 9mPi[—i{'m + I)] converges to a), 0,00 ', it follows 
that the second term in bracket ( |l05| ) tends to zero as -/V — > oo. Concerning the first term, we may 
write the inequality: 



N 



£ — ' to! 



m=0 



-/ I m + — 



< e 



N 1 

y — 



m=Q 



2 



(106) 



where the inequalities 



i/m 



< e, (to = 0, 1, 2, . . . , N) have been used. Next, by rewriting the 



Pollaczek polynomials Pe{—i(m + 1/2)] as 



Pi 



- i | m + - 



(107) 



and, substituting this expression in inequality (106), we obtain: 



N 1 

rn.\ 



E 

3=0 



TO+- 



(108) 



Next, we perform the limit for N — > oo. In view of the fact that X)j=o^^( TO + 1/2)^ is finite, and 
the series X)m=o( m + l/2) : '/ m ' converges, we can exchange the order of the sums and write: 



j=0 m=0 



(109) 



Finally, performing the limit for e — > 0, and recalling equality (105), statement (iii) is obtained, 
(iv) From definition (p9) it follows, for k% = fcp + 1, that: 



Eh 



> C. 



(110) 
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Statement (iv) (formula (|97|)) is proved if we can show that lim^y^^ k\ (e, N) = +00. Let us suppose 

that limjy^^ fci(e, N) is finite. Then there should exist a finite number K (independent of e and 

N) s uch that, for N tending to 00 and e tending to zero, k\(e,N) < K. Then, from inequality 
( |llfj| ) we have: 



C < 



fei(e.JV) 

E 

1=0 



K 

^E|« 

£=Q 



(e,JV) 



(111) 



But as N — > 00, e — > we have (recalling also statement (iii) formula (19 



< E |4 ' oo) 
^=0 



(112) 



which leads to a contradiction. Then statement (iv) follows. 

(v) Concerning statement (a), it follows obviously from definition (|98|) of M[ e ' N '. Finally, the first 



relationship in (|99|) is obvious; the second one follows from the asymptotic behavior of Pi[—i(m 
1/2)] at large I (for fixed m), i.e. formula (102). □ 



Remark 1. From statement (v) and formula ([)?]) it follows that the sum M£' N ' presents, for large 
values of N and small values of e, a plateau for k ~ fcrj. 



By truncating expansion (68) we may now introduce an approximation of the function F*(v) of 
the following type: 



F^ N \v) = £ a [ l' N) Mv), (v G R + ). 



(113) 



Approximation (v) is defined through the truncation number fcn(e, AT); the latter can be 

numerically determined by plotting the sum M^' N ^ versus k, and exploiting properties (a) and 
(b), proved in statement (v) of the previous lemma and the property stated in the remark above 
(see also Ref. [§). 

Now, we want to prove that the approximation F* N \v) converges asymptotically to F#(v) in 
the sense of the L 2 -norm, as N — > 00 and e — > 0. We can prove the following theorem. 

Theorem 5. The equality 



hm 

N—>oo 



F — F 

J- -Jf. J- 



(«) 



L 2 [0, + oo) 



(114) 



holds true. 

Proof. From the Parscval equality it follows that: 



F — F 



(e.JV) 



L 2 [0, + oo) 



E I- 

=fc +i 



(0,00) 



Eh 



(e,JV) (0,oo) 

Kin 



(115) 
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(0,oo) 



Since J2e=o 

0. It is convenient to rewrite the second term of the r.h.s. of (115) as follows. Let us define: 



= C and liirijv^^ fco(e, N) = +00, it follows that hin^^oo Efcfc +1 
£-►0 , , . -o 



(0,00) _ J af'°°^ if £ is even, 
I — ial 0,00 ' if I is odd, 



(116) 



[ (e,N) -r d ■ 

^(e,Af) _ I aj it i is even, 

I —iaf' N ^ if € is odd. 



(117) 



Notice that hf' 00 ^ and fif' N ^ are real, and Y^tLo 
Next, we introduce the following functions: 



(e,JV) ,(0,oo) 



ff (0>oc) (u)= g /l (0 ) oc) 1[ ^ +i[(?;) 



«=0 



^w = E4 e,iV) i[w[( 



(118) 
(119) 



e=o 



where 1^; is the characteristic function of the set E. From statements (i), (ii) and (hi) of the 
previous lemma (formulae (93), (|94|) and (Bq)) we obtain: 



p-\~OC 9 OO q 

/ (H^\v)) dv = £(/4 <°°>) = C 
' 1=0 



(120) 



-00, 



(121) 



^'^(u) — ► H {0 >°°Xv), (v G [0, +00)). 

JV^oo 
c^0 



(122) 



Hereafter, we assume, for the sake of simplicity and without loss of generality, that every term tif 

is different from zero. Next, let V(e, N) be the unique root of equation (H( e ' N \v)) dv = C. 

Let us indeed observe that (if («)) dv is a continuous non-decreasing function which is 
zero for V = 0, and +00 for V — ► +00. Furthermore, from statement (iv) of the previous lemma 
(formula (p7|)) we have limjy^^ F(e, AT) = +00. 

e-»0 

Then we can write: 



V(e,IV) 



(#(0,oc) ( 



V(e,iV) 



2 /-V(e,JV) f -i 

r)]dx-2J H^°°)(v) \H^' N \v) -/r(° ,0 °)(u)J dv. (123) 
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Next, we perform the limit for N ~* oo and e — > 0. Concerning the first term at the r.h.s. of (123) 
we have: 



lim 

N^oo Jvte,N) 



dv = 0. 



For what concerns the second term, we introduce the following function: 



B^ N \v) = { H(e ' N) ( v ) - H^(v) if < v < V(e, N) 



Then, we have by the use of the Schwarz inequality 



D 



(v) 



dv < AC, 



ifv> V(e,N). 



(N < oo,e > 0). 



Moreover, from (|l22j) we have: 



N^oo 



, v e [0, +oo). 



(124) 



(125) 



(126) 



(127) 



The family of functions {B^' N ^(v)} is bounded in L 2 [0, +oo), therefore it has a subsequence which 
is weakly convergent in L 2 [0, +oo). The limit of this subsequence is zero. In fact, let us observe that 
\B^ e ' N '(v)\ < 2C; then we consider the function B^ e,N ^ (v)<p(v), where cf> is an arbitrary element 
of the class of functions C^°(R + ). We then have \B^' N ^> (v)<f>(v)\ < 2C\(p(v)\, and this inequality 
does not depend on N an d e. In view of the Lebesgue dominated convergence theorem we can then 
write (see also limit (127)): 



lim sup 



N- 



B^ N \v)</>(v) dv 



= 0. 



(128) 



Since the set of functions C^°(R + ) is everywhere dense in L 2 [0, +oo), given an arbitrary function 



iji € i 2 [0,+oo) and an arbitrary number r\ > 0, there exists a function € C° 
II ^ ~~ 0fc II £2 [o +oo) < r l- Furthermore, through the Schwarz inequality we have: 



+ OC 



B^ N \v){Mv) - Hv)} 



dv < 

< 2VCrj. 

From fll28[) and (|129|) we can conclude that 



B 



(e,JV) 



(v) 



dv 



1/2 



\<l>k{v) - i>{v)\ dv 



lim sup 

]V->oc 



B^ N \v)ijj{v) dv 



0. 



) such that 

(129) 
(130) 



for any ip <E L 2 [Q, +oo). 

Next, by using the same type of arguments, we can state that if there is an arbitrary subsequence 
belonging to the family {B^ e ' N ^} that weakly converges in L 2 [0, +oo), then the weak limit of this 
subsequence is necessarily zero. Finally, from the uniqueness of the (weak) limit point, it follows 
that the whole family {bW*} converges weakly to zero in £ 2 [0, +oo). 
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We can thus write: 



+00 



lim / H {0 ^(v)B^ N \v)dv = 0, 
and from equality (|123[) we have 



lim 

— >oo In 



V(e,N) 



dv = 0. 



Since Y,iLn 



(e,JV) (0,oo) 



< 



jV(e,N) [ H (e,N)( v ^ _ ff(0,°°) („)] 2 dv > we have; 



fen 



lim } a 



(e,JV) (O,oo) 



= 0, 



and, in view of equality (115), the theorem is proved. □ 

We can then prove the following corollary. 
Corollary 1. The following equality holds true: 



lim 

N— >oo 
e->0 



k {e,N) 



Lf m) [0,+oo) 



Proof. The statement follows immediately from Theorem ^ by noting that: 



+00 



feo(e.JV) 



dv = 



(131) 



(132) 



(133) 



(134) 



L^ w \ [0, +00) feeing i/ie weighted L 2 -space with weight function w(v) = e v , (v G K + ), and the 
functions <l > £(v) are defined by formula (fzol). 



c/t' 



+00 k (e,N) 

F{v)-e v ' 2 4 e ' N) $t(v) 

k (e,N) 

F(v)-e v / 2 J2 4' N) Mv 



£=0 



(135) 



We can thus conclude that the jump function J^' b \v,-) — F(v) can be approximated by the 
truncated expansion 



fcoO.AT) 

£=0 



□ 



(136) 
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